Disordering Transitions in Vortex Matter: Peak Effect and Phase Diagram 



C. J. Olson, C. Reichhardt, R. T. Scalettar, and G. T. Zimanyi 

Department of Physics, University of California, Davis, California 95616. 

Niels Gr0nbech- Jensen 
Department of Applied Science, University of California, Davis, California 95616. 
NERSC, Lawrence Berkeley National Laboratory, Berkeley, California 94720. 

(February 1, 2008) 

Using numerical simulations of magnetically interacting vortices in disordered layered supercon- 
ductors we obtain the static vortex phase diagram as a function of magnetic field and temperature. 
For increasing field or temperature, we find a transition from ordered straight vortices to disordered 
decoupled vortices. This transition is associated with a peak effect in the critical current. For 
samples with increasing disorder strength the field at which the decoupling occurs decreases. Long 
range, nonlinear interactions in the c-axis are required to observe the effect. 
PACS numbers: 74.60. Ge, 74.60.Jg 



Vortex matter in superconductors exhibits a remark- 
ably rich variety of distinct phases due to the numerous 
competing interactions JT|. These phases can strongly 
affect the response of the system, such as the critical 
current J c and magnetization. In highly anisotropic su- 
perconductors, such as BSCCO, a pronounced fish-tail 
or peak effect, in which J c shows a sharp increase, is ob- 
served as a function of increasing field §. This peak 
can be interpreted as occurring when the increasing field 
weakens the interlayer coupling of vortex pancakes due 
to geometric constraints, and a transition occurs from 
weakly pinned 3D line vortices to decoupled 2D pancake 
vortices which can more easily adjust their positions to 
maximize the pinning The peak has also been pro- 
posed to arise from plasticity, proliferation of in-plane 
defects, dynamical effects, or matching effects 

There is mounting experimental evidence that the peak 
effect is associated with a sharp transition in the vor- 
tex lattice from an ordered state to a disordered state. 
In BSCCO, neutron scattering JjJ and muon lifetime 
Q experiments provide evidence that a transition from 
an ordered 3D vortex arrangement to a disordered or 
decoupled arrangement is associated with the peak ef- 
fect. Additional evidence from plasma-resonance j^] and 
magneto-optical studies point to the first-order na- 
ture of this transition. In YBCO a rapid increase in J c 
as a function of magnetic field is observed and is thought 
to indicate a transition from an ordered to a disordered 
state. History and memory effects found near this peak 
indicate that this transition is also first order ||, sug- 
gesting that the physics of the second peak is similar in 
BSCCO and YBCO. In recent muon-spin measurements 
in YBCO, the order-disorder transition has been inter- 
preted to be a 3D-2D transition of the vortex lattice [[ll| . 
As a function of temperature, in YBCO a peak effect 
can be observed very near or at the melting line fll2| , |i"3|] . 
Transformer measurements in this regime provide evi- 
dence of vortex cutting in the liquid state, suggesting 
that the breakdown of the 3D nature of vortex lines is 



relevant in this case as well. 

A key question is what type of sharp transition is re- 
sponsible for the observed peak in J c , and whether the 
mechanism of the peak effect is the same going through 
the melting line as through the second peak. In this pa- 
per we argue that it is a decoupling transition along the 
vortex line, in the c-direction, which is responsible. We 
present results from a simulation of magnetically inter- 
acting vortices in which we demonstrate a robust peak 
in J c as a function of magnetic field and produce an H-T 
phase diagram in good agreement with experiment. 

We consider a 3D layered superconducting model ma- 
terial containing an equal number of pancake vortices in 
each layer, interacting magnetically. This simulation dif- 
fers from those used previously by several other groups in 
that it treats the long-range interactions along the length 
of the vortex line exactly, and neglects the Josephson cou- 
pling, approximating highly anisotropic materials. Pre- 
vious models have treated the inter-plane interactions as 
simple (linear) elastic connections [[ll] pT[ . 

The overdamped equation of motion for vortex i is fj = 
- Ef=i ViU(p - , z tJ ■) + ff + f d + f T = Vi. The total 
number of pancakes per layer is N v , and pij and are 
the distance between vortex i and vortex j in cylindrical 
coordinates. We impose periodic boundary conditions in 
x and y directions The interaction energy between 
pancakes is |l9|,p0| 
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where R is the maximum in-plane distance, 
J p 00 d /9 'exp(pVA)/p', E 2 = /^dp'expCv^T^/A)/^, 
e = $q/(47tA) 2 , d — 0.005A is the interlayer spacing, 
and A is the London penetration depth. We model the 
pinning as N p short range attractive parabolic traps that 
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are randomly distributed in each layer. The pinning in- 



teraction is f- p = - T,k=i(fp/€p)( T i - r D e (fe> - l r » - 

|)/A), where the pin radius is £ p , the pinning force is 
fp, and /q = eo/A. Thermal fluctuations are represented 
as Gaussian noise of width fx, with < ff(t) >= and 
< fi{t)fj{t') >= 2r)k B TS l: jS(t - t'). The parameter 
s m is used to vary the coupling strength between planes, 
but except where noted it is set to s m = 1.0. To vary the 
applied magnetic field H , we fix the number of vortices 
in the system and change the system size (vortex den- 
sity n v ). The pin density remains fixed, but in all cases 
N p > N v . We consider systems of L = 8 to 32 layers 
containing from 1 to 80 vortex pancakes per layer. The 
system size ranges from 0.6A x 0.6A to 141A x 141A. 

The pancakes in our model can behave in two pos- 
sible ways. They may align along the c-axis into well- 
defined vortex lines ("3D"). Alternatively, the pancakes 
may break apart in the c-direction and move indepen- 
dently in each plane ("2D"). We can cross between the 
two types of behavior in the presence of disorder by 
varying the inter-plane coupling strength, s m . The vor- 
tices are 3D at high coupling, while 2D at low cou- 
pling. We find that due to the long range interac- 
tions and the nonlinearity, the transition between these 
two states is sharp and has many first-order character- 
istics, including strong hysteresis as well as superheat- 
ing and supercooling effects. We have studied this cou- 
pling strength crossover in detail in previous papers [^l| , 
and have shown that it is associated with a large change 
in the critical current / c , as indicated in Fig. 1(a). To 
quantify the alignment of the vortices between layers, 
we also plot the correlation function in the z-direction, 
C z = 1 - (9(a /2 - |r^ L - r^^+i | ) |r l:i - r^+i |2/a ), 
where ao is the vortex lattice constant. When C z = 1.0, 
the pancakes are aligned into 3D lines, whereas a low 
value of C' z indicates that the pancakes are decoupled. 

It is expected that changing the magnetic field in the 
system will cause a decoupling transition for the same 
reason that changing the coupling strength for inter- 
layer interactions does. This is because the vortices 
come closer together in plane as the field increases while 
the spacing between planes remains constant. Thus, at 
higher fields the relative coupling between planes be- 
comes weaker. Here, we show that the same decoupling 
transition previously observed as a function of coupling 
strength also occurs as a function of magnetic field. 

In Fig. l(b-d) we plot f c as a function of vortex density 
(magnetic field) n v , at a fixed temperature of fx = 0.0005 
for several different samples containing different numbers 
of vortices. A clear peak or fish-tail effect is present. In 
Fig. 1(c), a sample containing N v — 4 is shown. Here, for 
low fields (n v < 0.1), f c is high and the vortices are un- 
corrected in the z direction as well as disordered in plane. 
At these low fields the vortex-pin interactions dominate. 
f c drops rapidly with field in this regime due to the 
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FIG. 1. (a) 3D-2D transition as a function of inverse inter- 
layer coupling strength, l/s m , for a sample in which N v — 80, 
L = 16, n v = 0.35, n v = 1.0, and f p = 0.02/ *. Filled cir- 
cles: / c ; open squares: C z - (b-d) f c (filled circles) and C z 
(open squares) as a function of vortex density n v in samples 
with s m = 1.0 and: (b) N v = 1, L = 32, n p = 40.0, and 
f p = 0.2/ *. (c) N v = 4, L = 8, n p = 2.0, and f p = 0.04/ *; 
(d) N v = 80, L = 16, n p = 8.0, and f p = 0.1/ * . 

increasing in-plane vortex lattice stiffness, which causes 
the pinning to be less effective p^| . For intermediate 
fields, 0.1 < n v < 0.2, f c is low and the vortices form an 
ordered 3D structure as indicated by a near unity corre- 
lation, C z ~ 1. The 3D vortex lines in this regime are 
poorly pinned by the point pinning. At n v = 0.1 there 
is a sharp decoupling transition as is apparent in the 
abrupt drop of C' z . Simultaneously, f c increases rapidly 
to f c = 0.8/p. At higher fields f c gradually decreases due 
to the increasing in-plane interactions, which continue to 
stiffen the vortex lattice and thereby weaken the effec- 
tiveness of the pinning. We have observed similar peaks 
in the critical current for samples containing N v = 12, 
30, and 80 vortices. Data for the N v = 80 sample is 
shown in Fig. 1(d). The exact value of n v at which the 
peak falls is affected by the pinning strength and can be 
adjusted within a wide range, as indicated in Fig. 1. 

To show that the peak effect observed here arises from 
the onset of plasticity of the vortices between layers, and 
not from in-plane plasticity, we conduct simulations of a 
rigid vortex lattice, represented by a single vortex with 
periodic boundary conditions. In the N v = 1 simulation 
in-plane defects are not permitted. Thus, the vortex line 
can disorder along the c-axis only. The fact that we can 
observe the second peak in a sample containing only a 
rigid lattice, as in Fig. 1(b), indicates that the observed 
peak is not being generated by the proliferation of defects 
in the plane. This, along with our previous results on the 
existence of a transition as a function of c-axis coupling 
strength, provides evidence that it is a change in c-axis 
correlation that causes striking changes in f c . 

We can now probe the effect of temperature on the 
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FIG. 2. (a) f c versus n v for a sample with N v = 4, L — 8, 
and / p = 0.04/o at temperatures of fr = 0.0005 (circles), 
/t = 0.1 (squares), fx = 0.2 (diamonds), and fr = 0.3 (tri- 
angles), (b) C z corresponding to panel (a), (c) V x versus fr 
for the same sample, at a fixed drive of fd = 0.1/q , for fields 
n v = 0.061, 0.083, 0.108, 0.133, and 0.178. (d) C z correspond- 
ing to panel (c). Inset to (c): V x versus fr for a sample with 
N v = 80, n v = 2.0, L = 16, n p = 8.0, and f p = 0.1/,*. 

observed transition. As illustrated in Fig. 2(a-b), the field 
at which the 3D-2D transition occurs remains almost con- 
stant with T, in agreement with experiment ]23| ]. f c in 
the decoupled state drops as the temperature is raised, 
but the jump in f c persists. In Fig. 2(c) we show the 
peak effect in the average vortex velocity in the direc- 
tion of drive, V x =< v x >, as a function of temperature 
for fixed values of n v . Starting at fx = 0, we apply a 
constant driving force of fd = 0.1/q such that the vor- 
tices are moving and ordered. As fx is increased there 
is a transition to a 2D decoupled vortex arrangement as 
seen in the drop in C z . This coincides with a sharp drop 
in the vortex velocities V x since the decoupled vortices 
experience stronger effective pinning. As fx is further 
increased, the effectiveness of the pinning gradually de- 
creases and V x gradually increases back to the free flow 
level. The f c measurements at different temperatures 
and the resistivity measurements as a function of tem- 
perature exhibit the same features as experimental data 
taken near the peak effect regime. 

Next, we construct a phase diagram for our model in 
the H-T plane. Since we wish to stress that all qualita- 
tive features of the phase diagram result from the c-axis 
interactions, we consider a system containing a rigid vor- 
tex lattice. We first locate the thermal melting line for 
N v = 1 in the absence of pinning. This is shown as open 
circles on the phase diagram. We find a reentrance in 
the melting, shown in more detail in the top inset of Fig. 
3, which agrees with theoretical predictions for the melt- 
ing In the presence of pinning, we measure f c as a 
function of field (as in Fig. 1) at several different tem- 
peratures. We plot the 3D-2D transition field as solid 
squares. As was shown in Fig. 2, the location of the 
transition is insensitive to temperature at low 




FIG. 3. Phase diagram for a sample with N v = 1 and 
L — 32. Open circles: decoupling transition for f p = 0.0 
(clean melting line). Filled squares: decoupling transition for 
a sample with f p = 0.2/q. Triangles: decoupling for a sam- 
ple with f p = 1.6/q. Upper inset: Reentrance in the clean 
melting line is shown on a log- log plot. Lower inset: Phase 
diagram for a sample with iV„ = 4 and L = 8, with stronger 
pinning of f p — 0.04/q*. Circles: clean melting line. Squares: 
decoupling transition determined from f c vs. n v data as in 
Fig. 2(a). Plus signs: decoupling transition determined from 
V x vs. fr data as in Fig. 2(c). 

temperature. We find, however, that the temperature 
weakens the pinning enough that for /t > 1 a transition 
in f c is no longer observed. Instead f c is small every- 
where. The transition is still observable in C z , and data 
taken from C z is shown above fx = 1, indicating that the 
transition begins to drop in field at higher temperatures 
until it merges with the clean melting line. 

The field at which the second peak occurs is lowered 
in samples with stronger pinning, as indicated in Fig. 3 
where the second peak line for a sample with N v = 1 
and stronger pinning is plotted. This is consistent with 
experiments in which the peak shifts to lower fields in 
samples that have had their pinning strength increased 
artificially. 

We have also constructed a phase diagram for a sam- 
ple with N v — 4 and much stronger pinning, shown in 
the lower inset of Fig. 3. In the disordered system the 
3D-2D transition is relatively flat in fx for fx < 0.3 and 
then begins to turn down as the clean melting line is ap- 
proached. We also observe that the clean system shows a 
reentrant behavior for low fields due to the reduced vor- 
tex interactions. In the disordered sample this reentrant 
line is pushed up in field. We see reentrance of the second 
peak line in good agreement with recent measurements 
p3| . In still larger systems containing N v = 80 vortices, 
we have taken two slices across the phase diagram as a 
function of field and temperature, shown in Fig. 1(d) and 
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FIG. 4. Critical current as a function of magnetic field for a 
model containing springs between layers, (a) N v = 1, L = 32, 
fp = 0.2/ *. (b) N v = 4, L = 8, f p = 0.04/ *. (c) iV„ = 12, 
L = 16, ,f p = 0.1/q*. (d) iV„ = 80, L = 16, /„ = 0.1/ *. 

the inset of Fig. 2(c), respectively, and find decoupling 
transitions consistent with those shown above. 

To show that it is the long range interactions and the 
nonlinearity in the c-axis that cause the sharpness of the 
transition, we replace the long range interplane interac- 
tions by a model in which the pancakes are coupled by 
elastic nearest-neighbor springs. This model resembles 
one considered earlier |16j^]. The interlayer interaction 
is f it j = —s m (r iiL+ i - 2r^ L + r i>L _i), where s m can be 
used to vary the coupling strength. 

As shown in Fig. 4, for varying numbers of vortices, 
we hnd a smooth change in depinning force for the same 
parameters and over the entire range of n v considered in 
the long range interacting case. For N v = 1, the inter- 
layer interaction of the spring is not affected by varying 
n v so we vary s m instead. For N v = 4, 12, and 80, we 
set s m = 1 and vary n v directly. In each case we hnd a 
smooth decrease of f c with held. We do not observe any 
indications that a model containing springs and lacking 
long range interaction and nonlinearity in the z-direction 
can show a statistically significant peak in the critical 
current, as was recently suggested p"7| , p5| . 

In conclusion, using numerical simulations of magnet- 
ically interacting vortex pancakes in three dimensional 
layered superconductors, we have obtained static vortex 
phase diagrams as a function of magnetic held and tem- 
perature. We have demonstrated a unique relationship 
between c-axis vortex correlations and the critical current 
J c . Specifically, we have identified c-axis correlation tran- 
sitions to be responsible for the peaks in J c as well as the 
hsh-tail effect in magnetization measurements. Our sim- 
ulations of rigid in-plane vortex lattices, which also ex- 
hibit the peak and hsh-tail phenomena, support the claim 
that changes in c-axis (as opposed to in-plane) correla- 
tions are responsible for the observed anomalies. Finally, 
we have addressed recent suggestions that the peak effect 



in J c can be observed in a model system where inter- 
plane interactions are modeled by linear springs between 
pancakes. We only observe the critical current peak when 
long-range interactions are present, again suggesting that 
the peak and hsh-tail effect are directly related to the 
nonlinear and nonlocal nature of the inter-plane vortex 
interactions that lead to decoupling transitions. 
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